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Abstract 

We study the influence of the dimension of space on the thermodynamics of the classical and quantum self- 
gravitating gas. We consider Hamiltonian systems of self-gravitating particles described by the microcanonical 
ensemble and self-gravitating Brownian particles described by the canonical ensemble. We present a gallery of 
caloric curves in different dimensions of space and discuss the nature of phase transitions as a function of the 
dimension d. We also provide the general form of the Virial theorem in d dimensions and discuss the particularity 
of the dimension d = 4 for Hamiltonian systems and the dimension d — 2 for Brownian systems. 



1. Introduction 



Self-gravitating systems have a strange thermodynamics due to the attractive long-range nature of the grav- 
itational interaction [1,2,3,4]. This gives rise to inequivalence of statistical ensembles associated with regions of 
negative specific heats in the microcanonical ensemble and to phase transitions associated with gravitational col- 
lapse (see, e.g., [5]). The form of the caloric curve and the nature of the phase transitions crucially depend on the 
dimension of space d [6,7]. In this paper, we study the influence of the dimension of space on the thermodynamics 
of self-gravitating systems. We distinguish between Hamiltonian and Brownian systems [8]. A Hamiltonian system 
of self-gravitating particles evolves at fixed mass and energy and is described by the microcanonical ensemble [1]. 
In the mean-field approximation valid in a proper thermodynamic limit N — > +oo with fixed A = —ER/GM 2 
(where R is the size of the system), the distribution function maximizes, at statistical equilibrium, the Boltzmann 
entropy at fixed mass and energy. The relaxation towards statistical equilibrium is governed by the mean-field 
Landau-Poisson system which conserves mass and energy and monotonically increases the Boltzmann entropy 
(//-theorem). A system of self-gravitating Brownian particles [9] evolves at fixed mass and temperature and is 
described by the canonical ensemble. In the mean-field approximation valid in a proper thermodynamic limit 
N — > +oo with fixed rj = flGMm/R, the distribution function minimizes, at statistical equilibrium, the Boltz- 
mann free energy at fixed mass. The relaxation towards statistical equilibrium is governed by the mean-field 
Kramers-Poisson system which conserves mass and monotonically decreases the Boltzmann free energy. 

These two systems (Hamiltonian and Brownian) are defined in Sec. 2 starting from A-body equations. In Sec. 
3, we derive the general form of Virial theorem in d dimensions and we discuss the particularities of the dimension 
d = 4 for Hamiltonian systems and the dimension d = 2 for Brownian systems. In the first case, we show that 
the system evaporates for E > and collapses for E < 0. In the second case, the system evaporates for T > T c 
and collapses for T < T c where T c = (N - l)Gm 2 /4k B is a critical temperature. We also derive a generalization 
of the Einstein relation including self-gravity in d — 2. Finally, in Sec. 4, we present a gallery of caloric curves in 
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various dimensions of space for self-gravitating classical and quantum particles (fermions) and discuss the nature 
of phase transitions as a function of d. 



2. Self-gravitating Hamiltonian and Brownian systems 



In a space of dimension d, we consider a system of N particles with mass m a in gravitational interaction whose 
dynamics is described by the equations of motion 
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Here, the greak letters refer to the particles and the latin letters to the coordinates of space. When the last two 
terms are set equal to zero, we recover the usual Hamiltonian model of self-gravitating systems. However, we 
consider here a more general situation where the particles are subject, in addition to self-gravity, to a friction 
force and a stochastic force. These terms can mimick the influence of a thermal bath of non-gravitational origin. 
We thus obtain a model of self-gravitating Brownian particles [9] extending the usual Brownian model introduced 
by Einstein and Smoluchowski. We consider the case of a multicomponent system [10]. Here, £ a is the friction 
coefficient, D a is the diffusion coefficient and R™(t) is a white noise acting independently on the particles and 
satisfying (Rf(t)) = and {R?(t)Rj(t')) = 8ij8 a [j8(t — t'). The diffusion coefficient and the friction force satisfy 
the Einstein relation D a = ^ a ksT /m a where T is the temperature of the bath (see below). Since D ~ T, the 
temperature measures the strength of the stochastic force. 

The evolution of the JV-body distribution function Pjv(ri,vi, rjv, vjv, t) is governed by the iV-body Fokker- 
Planck equation [8]: 
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where F a — — V a f/(ri, ...,rjv) is the force acting on the a-th particle (17 is the potential of interaction). If we 
enclose the system within a box and replace the bare gravitational potential by a soften potential that is regularized 
at short distances (so that a statistical equilibrium state exists), we find that the stationary solution of Eq. (2) is 
the canonical distribution 



Pjv(ri, vi, ...,rjv,vjv) = 
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where H — X] a m <*"2 L + t^( r i, r iv) is the Hamiltonian, provided that the diffusion and friction coefficients 
are related to each other by the Einstein relation. When £ = D = 0, Eq. (2) becomes the Liouville equation 
appropriate to Hamiltonian systems. In that case, Eq. (3) is not valid anymore. Since energy is now conserved 
during the evolution, the system is expected to reach, for t — » +oo, the microcanonical distribution 
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expressing the equiprobability of accessible microstates (with the right value of energy). 

In the mean-field approximation valid in a proper thermodynamic limit with N — ► +oo [8,10], the evolution 
of the distribution function / a (r, v, t) of each species of the self-gravitating Brownian gas is governed by the 
multicomponent Kramers-Poisson (KP) system 
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where / = f a is the total distribution function. The multicomponent KP system conserves the mass M a — 
J fadrdv of each species and monotonically decreases the free energy F — E — TS where 
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is the total energy and 
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is the Boltzmann entropy. The linearly dynamically stable stationary solution of the KP system is the mean-field 
Maxwell-Boltzmann distribution 

/ (r,v) = A a e-^^™ (9) 

that is a local minimum of free energy at fixed mass (canonical description). In the strong friction limit £ — > +00, 
the velocity distribution of the particles is close to the Maxwellian and the evolution of the spatial distribution 
is governed by the multicomponent Smoluchowski-Poisson (SP) system [8,10]. Alternatively, for D = £ = 
(Hamiltonian systems), the evolution of the system for N — ► +00 is governed by the Vlasov equation [8]: 

The Vlasov equation conserves mass and energy but also the entropy (and more generally an infinite class of 
functionals called the Casimirs). Furthermore, it admits an infinite number of stationary states (not necessarily 
Boltzmannian) that can be selected by the dynamics as a result of an incomplete collisionless violent relaxation 
on the coarse-grained scale [11,12]. The statistical equilibrium state in the microcanonical ensemble is selected 
by finite N effects (granularities) accounting for correlations between particles due to close encounters. The 
collisional relaxation of stellar systems is usually described by the mean-field Landau-Poisson system [8] . In d = 3 
the multicomponent Landau equation can be written 
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where u = v — v' is the relative velocity of the particles involved in an encounter and In A = / + °° dk/k is the 
gravitational Coulomb factor. We have set /„ = / (r,v',£) assuming that the encounters can be treated as local. 
The Landau-Poisson system monotonically increases the entropy while conserving the mass M a of each species and 
the total energy E. The linearly dynamically stable stationary state is the mean-field Maxwell-Boltzmann distri- 
bution (9) that is a local maximum of Boltzmann entropy at fixed mass and energy (microcanonical description). 
Therefore, the statistical equilibrium states of Hamiltonian and Brownian systems have the same distribution 
profiles but their stability may differ in case of ensemble inequivalence. This occurs in particular when the caloric 
curve presents turning points as discussed in Sec. 4. 



3. Virial theorem in d dimensions 



3.1. The general expression 

Let us establish the Virial theorem associated with the stochastic equations (1). For convenience, we shall 
assume that the friction coefficient £ is the same for all the particles but this assumption can be relaxed easily. 
The moment of inertia tensor is defined by 

Iij = y^ y m ct xfx'j. (13) 

We introduce the kinetic energy tensor 

K i j = \52m a x?x?, (14) 
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and the potential energy tensor 

Wij = G^m a m Pl -—- F = - 2 G^m a m p ^— , (15) 

where the second equality results from simple algebraic manipulations obtained by exchanging the dummy vari- 
ables a and (3. Taking the second derivative of Eq. (13), using the equation of motion (1), and averaging over the 
noise and on statistical realizations, we get 

+ \t{ia) = 2(Kij) + (Wij) - i j (PikXj + P jkXi ) dS k . (16) 

For the sake of generality, we have accounted for the presence of a box enclosing the system and included the 
Virial of the pressure force Fi = PikASk on the boundary of the box: (Iij)box = 2~2box2~2a(Pi" x 7 + Fj" x ?) = 
2~2box EaK 1 ? + Pjk x i) ASk- The scalar Virial theorem is obtained by contracting the indices yielding 
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where 
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are the moment of inertia and the kinetic energy. On the other hand, 

2 ^ \r f 3-r a \ d 2 
3.2. The trace of the potential energy tensor 
For d 7^ 2, we find that 

Wu = (d-2)W (20) 
where W is the potential energy 

W = — T mamp (21) 

In that case, the scalar Virial theorem reads 

+ = 2(K) +(d- 2){W) - j P ik XidS k . (22) 

For d = 2, we have the simple result 

Wu = ~\g ^2 m a m/3. (23) 

For equal mass particles, 

Wu = -\gN{N -l)m 2 , (24) 
which reduces to 

Wu^-\CM 2 , (25) 

for N 3> 1. We also note that X! a ^/3 m a vrij) — M 2 — m 2 a . For typical mass distributions, the first term is of 
order ~ N 2 m 2 and the second of order Nm 2 (where m is a typical mass). Therefore, in the thermodynamic limit 
N — > +oo, we get Eq. (25) even for a multicomponent system [10]. 
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3.3. The equilibrium state 



At equilibrium, the Virial theorem (17) reduces to 

2{K) + (Wu) = j P tk x l dS k . (26) 

For an unbounded domain (P = 0), we get 2(K) + (Wu) = (more precisely, 2{K) = -(d - 2)(W) for d 2 
and (K) = GM 2 /4 for d — 2 in the approximation (25)). If the system is at statistical equilibrium, then (K) = 
^NUbT and Pij = pb~ij with p = J] s p s k B T ' jm s (where s labels the species of particles) [10]. If pt is constant on 
the boundary of the box, the pressure term becomes 

Pb v-dS=p b jv-dS=p b j V • r dV = dp b V. (27) 

More generally, we define P = § pr ■ dS. In that case, the equilibrium Virial theorem becomes 

dNk B T + (Wu) = dPV. (28) 
For an ideal gas without self-gravity (W = 0), we recover the perfect gas law 

PV = Nk B T. (29) 
Alternatively, for a self-gravitating gas in two dimensions, using Eq. (23), we get the exact equation of state 

PV = Nk B (T-T c ), (30) 
with the exact critical temperature 

B c ~ In • ( } 

For equal mass particles, we get 

k B T c = ^-(N-l), (32) 
and in the N — * +oo limit 

kBTc= ^N~- (33) 

This relation is also valid for a multicomponent system in the thermodynamic limit (see remark after Eq. (25)) 
[10]. The equation of state (30) has been obtained by Lynden-Bell & Katz [13], Padmanabhan [1] and Chavanis 
& Sire [14] using different methods. 

3.4. The critical dimension d = 4 for Hamiltonian systems 

For Hamiltonian systems (£) = £ = 0), the total energy E = K + W is conserved. Thus, in an unbounded 
domain (P = 0), the Virial theorem (17) can be written for d ^ 2: 

i I = 2K + (d- 2)W = (4 - d)K + (d - 2)E = 2E+{d- 4)W. (34) 

This is the extension of the usual Virial theorem in d dimensions (this equation is exact without averages). We 
note that the dimension d = 4 is critical as was also noticed in [7] using different arguments. In that case, I — 4E 
which yields after integration I — 2Et 2 + C\t + C2. For E > 0, I — > +00 indicating that the system evaporates 
(there can be equilibrium states if the system is confined within a box). For E < 0, I goes to zero in a finite time, 
indicating that the system forms a Dirac peak in a finite time. More generally, for d > 4, since (d — 4)W < 0, wc 
have I < 2Et 2 + Cit + C2 so that the system forms a Dirac peak in a finite time if E < (this result remains valid 
for box-confined systems). Therefore, self-gravitating systems with E < are not stable in a space of dimension 
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d > 4. The study in [7] indicates that this observation remains true even if quantum effects (Pauli exclusion 
principle for fermions) are taken into account. This is a striking result because quantum mechanics stabilizes 
matter against gravitational collapse in d < 4 [5]. For 2 < d < 4, since (d — 4)W > we conclude, according to 
Eq. (34), that the system evaporates if E > (there can be equilibrium states if the system is confined within 
a box). An equilibrium is possible, but not compulsory, if E < (the case of statistical equilibrium is discussed 
in Sec. 4). Finally, for d < 2, since W > 0, the energy is necessarily positive (E > 0). Since (d — 4)W < 0, an 
equilibrium state is possible. In d = 2, the Virial theorem becomes with the approximation (25): 

Since K > 0, an equilibrium state is possible. 

3.5. The critical dimension d = 2 for Brownian systems 

We consider a self-gravitating Brownian gas in the strong friction limit £ — > +oo. To leading order in l/£, the 
iV-body distribution is given by 

P N {ri,vi,...,T N ,v N ,t) = e- l3 ^ ma ^$ N (r 1 ,...,r N ,t) + 0(C 1 ), (36) 

as can be deduced from Eq. (2) by canceling the term in brackets. From this expression, we find that Pij = pSij 
with p = J2 s Psk B T/m 3 and (K l3 ) = ±iVfc B T5 l3 even for the out-of-equilibrium problem. From Eq. (16), we 
obtain the overdamped Virial theorem for a self-gravitating Brownian gas 

= ( w a) + Nk B TS tj - i j p(S ik Xj + S jk Xi) dS k . (37) 

We can obtain this result in a different manner. In the strong friction limit £ — » +oo, the inertial term in Eq. (1) 
can be neglected so that the stochastic equations of motion reduce to 

xf = n a m a ^fi'f } + V^R?(t), (38) 

where D' a — k B T)i a is the diffusion coefficient in position space and fi a = (^m a ) _1 is the mobility. Taking the 
derivative of the tensor of inertia (13) and using Eq. (38), we get 

hj = jWij + J2 m«v / 2£Uz?fl"W + x?R?{t)]. (39) 

' a 

Averaging over the noise using (xfRf) = y/2D' a 8ij and on statistical realizations, we recover Eq. (37). The scalar 
Virial theorem obtained by contracting the indices reads 

if</> = dNk B T + (Wu) - dPV. (40) 

In particular, in d — 2, using Eq. (23) we find that 

^(i) = 2Nk B (T-T c )-2PV, (41) 

where T c is defined in Eq. (31). For T > T c , an equilibrium state is possible in a box, while the system evaporates 
(/ — » +oo) in an unbounded domain (P = 0). For T < T c , the system creates a Dirac peak in a finite time (/ = at 
t = t end ). These dynamical evolutions are studied precisely in [6,15] by solving the Smoluchowski-Poisson system. 

3.6. The generalized Einstein relation including self-gravity in d — 2 

In an unbounded domain (P = 0), the Virial relation (41) reduces to 

1 -i(I) = 2Nk B {T-T c ). (42) 
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Defining the spatial dispersion of the particles by 
, 2, _ (Em a rj} _ (J) 

{ ' J2m a M ' {i6 > 
we can rewrite the foregoing relation in the form 

Integrating this relation, we obtain 

<r 2 >(i) = ^f^(l-T c /T)t + <r 2 >o. (45) 
This relation suggests to introducing an effective diffusion coefficient 

D eff (T) = ?^f(l-T c /T), (46) 
so that 

{r 2 )(t)=4D eff (T)t + (r 2 ) . (47) 

For T 3> T c when gravitational effects become negligible with respect to thermal motion, the diffusion coefficient 
is given by the celebrated Einstein formula D — Nk B T/£M. However, Eq. (46) shows that the diffusion is less and 
less effective as temperature decreases and gravitational effects come into play. Therefore, relation (4-6) provides 
a generalization of the Einstein relation to the case of self- gravitating Brownian particles in d = 2. For T > T c , 
we have a diffusive motion (evaporation) with an effective diffusion coefficient depending linearly on the distance 
(T — T c ) to the critical temperature. For T = T c , the effective diffusion constant vanishes D(T C ) = so that the 
moment of inertia is conserved. Finally, for T < T c the effective diffusion coefficient is negative, implying finite 
time collapse. In particular, (r 2 ) = for t en a = (r 2 )o/4jD(T)| (where (r 2 )o is calculated from the center of mass). 
This corresponds to the formation of a Dirac peak p(r) = MS(r) at r = containing the whole mass [15]. 



4. Thermodynamics of self-gravitating classical and quantum particles 

In this section we present a gallery of caloric curves for self-gravitating classical and quantum particles (fermions) 
at statistical equilibrium in different dimensions of space. From the series of equilibria, we can directly infer the 
thermodynamical stability limits of the system in microcanonical and canonical ensembles by using the turning 
point argument of Poincare. We also directly see the critical parameters (energy and temperature) below which 
the system collapses. The microcanonical ensemble (MCE) is appropriate to isolated Hamiltonian systems (in 
which case the control parameter is the energy E) and the stable states are maxima of entropy S at fixed mass M 
and energy E. The canonical ensemble (CE) is appropriate to dissipative systems like self-gravitating Brownian 
particles or self-gravitating gaseous systems in contact with a thermal bath of non-gravitational origin (in which 
case the control parameter is the temperature T) and the stable states are minima of free energy F = E — TS 
at fixed mass M. The series of equilibria /3(E) (representing critical points of the thermodynamical potentials) 
are the same for the two systems but their stability (related to the sign of the second order variations of the 
thermodynamical potentials) can be different. A turning point of energy is associated with a loss of microcanonical 
stability (for Hamiltonian systems) and a turning point of temperature is associated with a loss of canonical 
stability (for Brownian systems). Stability can be regained at the next turning point if the curve turns in the 
reversed direction (as is the case for fermions). We refer to Katz [16], Padmanabhan [1] and Chavanis [17,18,5] 
for more precise statements of these results and illustrations. In the following, the caloric curves are given for 
box-confined systems. This is necessary to have existence of statistical equilibrium states in d > 2 (otherwise the 
system evaporates). We also briefly discuss the case of unbounded systems in d = 1 and d — 2. 

4.1. Classical particles 

In d — 1, there exists statistical equilibrium states for any accessible value of energy E > in MCE and any value 
of temperature T in CE (see Fig. 1-a). Since the series of equilibria is monotonic, they are stable (global maxima 
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Figure 1. Series of equilibria for classical isothermal spheres in d = 1, d = 2 and d — 3 dimensions. 




Figure 2. Scries of equilibria for classical isothermal spheres with 2 < d < 10. 

of S in MCE and global minima of F in CE) and the ensembles are equivalent. There also exists stable statistical 
equilibrium states in an unbounded domain. The density profile is known analytically [19,6,15]. According to the 
Virial theorem, we have 2K = W and K = ^NksT so that the caloric curve is simply E — ^NksT [15]. 

In d = 2, there are statistical equilibrium states for any value of energy E in MCE and for T > T c — GMm/AkB 
in CE (see Fig. 1-a). Their density profile is known analytically (see, e.g., [6]). Since the curve (3(E) is monotonic, 
they are stable (global maxima of S in MCE and global minima of F in CE) and the ensembles are equivalent. For 
T < T c the self-gravitating Brownian gas (canonical description) undergoes an isothermal collapse and ultimately 
forms a Dirac peak [6]. In an unbounded domain, there are stable statistical equilibrium states in the microcanon- 
ical ensemble (Hamiltonian systems) for any value of the energy [20]. According to the Virial theorem, we have 
K — GM 2 /4 and K = NksT so that these equilibrium states have the same temperature T c = GMm/AkB, inde- 
pendently on their energy E. By contrast, in the canonical ensemble (Brownian systems), there is no equilibrium 
state in an unbounded domain as discussed in [15] and in Sec. 3.6. 

In d = 3, the series of equilibria forms a spiral (see Fig. 1-b). Since the series of equilibria presents turning 
points, the ensembles are not equivalent. There exists stable statistical equilibrium states in MCE for E > E c — 
-0.335GM 2 /R (Antonov energy); they have a density contrast R < 709. There exists stable statistical equilibrium 
states in CE for T > T c = GMm/2.52kBR (Emden temperature); they have a density contrast R < 32.1. They 
corresponds to long-lived metastable states (local maxima of S in MCE and local minima of F in CE) whose 
lifetime ~ increases exponentially with the number of particles [21]. The states with density contrast 32.1 < 
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R < 709 are stable in MCE and unstable in CE; they have negative specific heats. For E < E c a self-gravitating 
Hamiltonian system (e.g. a globular cluster) undergoes a gravothermal catastrophe [22] . There is first a self-similar 
collapse leading to a finite time singularity (the central density becomes infinite in a finite time) and ultimately 
the system forms a binary star surrounded by a hot halo [23]. This structure has an infinite entropy at fixed 
energy (see Appendix A of [6]). For T < T c a self-gravitating Brownian system undergoes an isothermal collapse 
[9]. There is first a self-similar collapse leading to a finite time singularity [6] and ultimately the system forms 
a Dirac peak in the post-collapse regime [24]. This structure has an infinite free energy (see Appendix B of [6]). 
There is no equilibrium state (no maximum of entropy or minimum of free energy) in an unbounded domain: the 
system evaporates or collapses [6,15]. 

In d > 3, the phenomenology is the same as in d — 3 (see Fig. 2) but we note, for curiosity, that the spiral 
disappears for cZ > 10 [6]. On the other hand for d > 4 and E < E c , the system forms a Dirac peak in the MCE 
instead of a "binary star + hot halo" structure [6] (see Appendix A of [6] ) . 

4.2. Self- gravitating fermions 

We now present the series of equilibria for the self-gravitating Fermi gas (described by the Fermi-Dirac entropy) 
in different dimensions of space and for different values of the degeneracy parameter The degeneracy parameter 
can be viewed as a normalized Planck constant [i ~ h~ d , as an effective inverse small-scale cut-off \jl ~ 1/e or as 
the system size /x ~ ^ d ( 4 ~ d )/ 2 (the classical case is recovered for h — > 0, e — > 0, R — > +oo or fj, — > +co). We refer 
to [5,7,25] for details. 

In d = 1, there exists statistical equilibrium states for any accessible value of energy E > E g (where E g is the 
ground state) in MCE and any value of temperature T in CE (see Fig. 3-a). They are stable (global maxima of 
S in MCE and global minima of F in CE). 

In d = 2, there exists statistical equilibrium states for any accessible value of energy E > E g in MCE and any 
value of temperature T in CE (see Fig. 3-b). They are stable (global maxima of S in MCE and global minima of 
F in CE). For /i — > +oo, E g — > — oo and we recover the classical caloric curve displaying the critical temperature 
T c . Below T c , the classical Brownian gas collapses and creates a Dirac peak ("black hole") [6]. When quantum 
mechanics is accounted for, the "black hole" is replaced by a "fermion ball" (or white dwarf star) consisting of a 
dense degenerate nucleus surrounded by a dilute atmosphere [7] . 

In d — 3, the nature of phase transitions depends on the value of the degeneracy parameter [5]. For high values 
of n (low small-scale cut-off), the series of equilibria has a Z-shape structure resembling a dinosaur's neck (see 
Fig. 4-a). Starting from a gaseous configuration (upper branch) and decreasing the energy, the system first passes 
from a stable gaseous phase (global maximum of S) for E > E t to a metastable gaseous phase (local maximum 
of S) for E < E t . The microcanonical first order phase transition at E t is avoided due to the long lifetime 
(~ e^) of the metastable gaseous states. At the critical energy E c , the gaseous metastable phase disappears 
and the system collapses (gravothermal catastrophe). However, the collapse is stopped by quantum mechanics 
(Pauli exclusion principle) when the core of the system becomes degenerate. Therefore, the system ends up in the 
condensed phase (lower branch). If we now increase the energy the system first passes from a stable condensed 
phase (global maximum of S) for E < Et to a metastable condensed phase (local maximum of S) for E > E t . 
Again, the microcanonical first order phase transition at E t is avoided due to the long lifetime of the metastable 
condensed states. At the critical energy E*, the condensed metastable phase disappears and the system explodes 
and returns to the gaseous phase. We can thus describe an hysteretic cycle in the MCE. For low values of n (large 
small-scale cut-off), the series of equilibria /3(E) has a iV-shape structure (see Fig. 4-b). Since there is no turning 
point of energy, all the states are stable in MCE (global maxima of S), for Hamiltonian systems, and there is no 
special phase transition for this value of degeneracy parameter, just a condensation (clustering) of the system as E 
decreases. Negative specific heats are possible in MCE. By contrast, since the curve E(0) has a Z-shape structure, 
there is a phase transition in the CE, for Brownian systems. Starting from a gaseous configuration (left branch) 
and decreasing the temperature, the system first passes from a gaseous stable phase (global minimum of F) for 
T > Tt to a gaseous metastable phase (local minimum of F) for T <Tt. The canonical first order phase transition 
at Tt is avoided due to the long lifetime (~ e") of the metastable gaseous states. At the critical temperature T c , 
the gaseous metastable phase disappears and the system collapses (isothermal collapse). However, the collapse 
is stopped by quantum mechanics (Pauli exclusion principle) when the core of the system becomes degenerate. 
Therefore, the system ends up in the condensed phase (right branch). We get a "fermion ball" instead of a "Dirac 
peak" for classical particles. If we now increase the temperature the system first passes from a stable condensed 
phase (global minimum of F) for T < T t to a metastable condensed phase (local minimum of F) for T > T t . Again, 
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Figure 3. Series of equilibria for an isothermal gas of self-gravitating fermions in d — 1 and d — 2 dimensions for different values of 
the degeneracy parameter (various system sizes). 




Figure 4. Series of equilibria for an isothermal gas of sclf-gravitating fermions in d — 3 dimensions for large (left, Z-shape) and 
small (right, TV-shape) values of the degeneracy parameter. 



the canonical first order phase transition at Tt is avoided due to the long lifetime of the metastable condensed 
states. At the critical temperature T*, the condensed metastable phase disappears and the system explodes and 
returns to the gaseous phase. We can thus describe an hysteretic cycle in the CE [26]. We note that when quantum 
effects are taken into account, there exists an equilibrium state for any value of accessible energy E > E g in MCE 
and any value of temperature in CE [5]. 

In d > 4, the situation is different [7]. In that case, there is no equilibrium state for sufficiently small energies 
and temperatures (see Fig. 5). This is similar to the Antonov instability for classical particles in d — 3 but this 
now occurs for fermions. Therefore, quantum mechanics cannot stabilize matter against gravitational collapse in 
d > 4, contrary to what happens in d = 3. This is consistent with our previous observation that classical white 
dwarf stars (the T = limit of the self-gravitating Fermi gas) would be unstable in a space of dimension d > 4 
[14]. This is also similar to a result found by Ehrenfest [27] at the molecular level (in Bohr's model). The case of 
relativistic white dwarf stars in d dimensions is treated in [28]. 



10 



e 



18 




D=4 


Incomplete | 
polytrope i 


H=23 


15 






I 




12 






// 


/ 


9 






H=20 / / 




6 








| H=25 


3 






H=30 ; 


I 


n 











-0.7 -0.5 -0.3 -0.1 0.1 

A=-ER D 2 /GM 2 

Figure 5. Series of equilibria for ail isothermal gas of self-gravitating fermions in d — 4 dimensions for different values of the 
degeneracy parameter. 
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